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ABSTRACT. 

We present two rigorous results on the Sherrington-Kirkpatrick mean field model for spin 
glasses, proven by elementary methods, based on properties of fiuctuations, with respect to 
the external quenched noise, of the thermodynamic variables and order parameters. The 
first result gives the uniform convergence of the quenched average of the free energy in 
the thermodynamic limit to its annealed approximation, in the high temperature regime, 
including the assumed critical point (/? = 1 in our notations). The second result shows 
that the free energy can be expressed through a functional order parameter, of the type 
introduced by Parisi in the frame of the replica symmetry breaking method. The functional 
order parameter is implicitely given in terms of fiuctuations of thermodynamic variables. 



I Research supported in part by MURST (Italian Minister of University and Scientific 
and Technological Research) and INFN (Italian National Institute for Nuclear Physics). 
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1. INTRODUCTION. 

The Sherrington-Kirkpatrick mean field model [1] for spin glasses is defined through the 
Hamiltonian 



where ai, . . . , aN are Ising spins, with values ±1, coupled in a highly disordered way 
through the quenched variables Jij, assumed to be independent random variables with 
identical unit Gaussian distribution (for the sake of simplicity). The sum runs over all the 
A'"(A'" — l)/2 couples of spins. 

Due to its relevant physical interest as a prototype for disordered complex systems, 
this model has been the subject of intensive investigation in the last years [2]. By now, 
its general structure is well understood, at least from a qualitative point of view, and 
ingeniously described by the Parisi solution [3], originally obtained through the replica 
symmetry breaking method, and further confirmed by the cavity method [4] . According to 
the accepted view, the model is trivial in the high temperature regime {/3 < 1), where the 
thermodynamic variables coincide with their annealed approximation, in the thermody- 
namic limit — )■ oo. On the other hand, at low temperatures (/3 > 1), the model exhibits 
an extremely rich structure of approximate equilibrium states, separated by high barriers, 
and organized in a hierarchical system. As a consequence, the thermodynamic variables 
are expressed through a functional order parameter, introduced by Parisi [5], replacing 
the numerical order parameter of the approximate Sherrington-Kirkpatrick solution [1], 
obtained without replica symmetry breaking and violating positivity of the entropy. 

Rigorous results, by Aizenman, Lebowitz and Ruelle [6] , give a very detailed descrip- 
tion of the fluctuations of thermodynamic variables and their limiting behavior, as A" — > oo, 
in the high temperature regime (0 < < 1) (see also [7]), and interesting estimates on the 
free energy at low temperature. 

On the other hand, Pastur and Scherbina [8] have given a rigorous proof of the damping 
of the mean square fluctuations of the free energy, in the thermodynamic limit, at all 
fixed temperatures. They have also shown that the approximate Sherrington-Kirkpatrick 
solution, surely defective at high /3, follows necessarily from the assumption that a suitably 
defined order parameter, connected with the induced magnetization in an external random 
field, has vanishing mean square fluctuations in the thermodynamic limit (see also [9]). 
This result strongly supports Parisi theory. 

The study of fluctuations is clearly of great importance for this model. In this paper, 
we give a rigorous proof of two results, based on elementary properties of fluctuations with 
respect to the external quenched noise J. The first result gives the uniform convergence 
of the quenched average of the free energy to its annealed approximation, in the high 
temperature regime, including the point P = 1, which is the critical point, according to the 
accepted picture. The second result is based on the evaluation of the corrections necessary 
to go from the quenched average to the annealed average, where the variables J participate 
to the thermodynamic equilibrium. These corrections are expressed through fiuctuations 
of appropriate thermodynamic variables. A simple consequence is the emergency of the 
functional order parameter and the antiparabolic martingale equation of Parisi theory [3] . 
Therefore, we show that the free energy can be expressed through the functional order 




(1.1) 
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parameter, exactly as in Parisi theory. For further developments along this line we refer 
to [10]. 

The paper is organized as follows. In Section 2, we give a detailed description of 
the model and some elementary properties of the thermodynamic variables. The results 
on the high temperature behavior are given in Section 3, where we exploit only simple 
convexity properties and positivity of fluctuations. In Section 4 we introduce the marginal 
martingale method, equivalent to the cavity method. It is a very powerful technique for 
the study of the thermodynamic limit in mean field models. 

The functional order parameter and the antiparabolic martingale equation are intro- 
duced in Section 5. They allow a very simple expression of the thermodynamic variables 
of the system. The functional order parameter can be easily written in terms of fluctua- 
tions, which take into account the corrections necessary to go from quenched averages to 
annealed averages. 

Finally, Section 6 is devoted to conclusions and outlook for future developments. 
2. THE MODEL. ELEMENTARY PROPERTIES. 

We consider N sites, i = 1, . . . , N. In the thermodynamic limit we will let N ^ oo. The 
configurations of the system are given by 

a: {l,2,...,iV}9z^a, eZ2 = {-l,l}, (2.1) 

where cr^ are Ising spins. For each of the A^(A^ — 1) couples of sites i 7^ j, we 

introduce independent random variables Jij = Jji, i ^ j, identically distributed , called 
quenched variables. The ex's are mesoscopic random variables subject to thermodynamic 
equilibrium. The J's do not participate to thermodynamic equilibrium, but act as a kind 
of random environment on the cr's. For the sake of simplicity, we assume that the J's have 
unit Gaussian distribution with 

£;(J,,)=0, i?(4) = l, (2.2) 

where E denotes averages with respect to the J variables. 
The Hamiltonian of the model is given by 

ifiv(c7, J) = - — ^= JijOiaj, (2.3) 

where the sum extends over all the A^(A^ — 1) distinct couples. The square root is necessary 
in order to assure the right thermodynamic behavior for the extensive variables, as it will 
be clear in the following. We write it in the form ^/N — 1, equivalent to the customary 
one given by y/N. Introducing the inverse temperature /3 (in proper units), we can write 
the Boltzmannfaktor as exp{—(3H]sf{a, J)), and define the partition function Z]sf{f3, J) and 
the free energy F^ili, J) through 

Zn{(5,J)^ exp(-/3i?Ar(cT,J))=exp(-/3FA,(/3,J)). (2.4) 

CTl...CTiV 
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The associated Boltzmann state 0Jn,i3,j is given by 

iON,p,j{A) = ZNW,J)-' Yl AeM-PHN{a,J)), (2.5) 

<Ti...ajv 

where ^ is a generic function of the cr's. Its density is 

PN{a; 13, J) = Z^iP, eM-PH^ia, J)). (2.6) 
The internal energy Un{/3,J) is 

UN{P,J)=ujN,f5AHN{a,J)) = -dplogZN{f3,J) = dp{pFN{P,J)). (2.7) 
The entropy Sn{/3,J) is 

Sn{/3,J) = - Y1 (PivlogPiv)(a;/3,J), (2.8) 

CTl...CTiV 

with the usual bounds 

0< Sn{I3, J) < Nlog2. (2.9) 
These variables are connected through the second principle of thermodynamics 

Fn{P, J) = Un{P, J) - r^SNiP, J). (2.10) 

All thermodynamic variables are random with respect to the external noise J. However, 
we can also introduce the free energy per site fN{P, J), and correspondently the internal 
energy per site un{/3, J), and the entropy per site sn{P, J), such that 

/nW, J) = N-'FnW, J), un{P, J) = N-'Un{I3, J), sn{P, J) = N-^SnW, J), 
< snW, J) < log 2, /jv(A J) = unW, J) - r^SNiP, J). (2.11) 

On physical grounds, it is expected that these variables have a good thermodynamic be- 
havior in the limit when A?" — > oo, in the sense that they should converge almost surely, 
with respect to the external J noise, i.e. with the possible exclusion of a set of J samples 
of zero measure in the J probability space, 

fN{l3,J)^ fNiP), Un{P,J)^Un{I3), SNiP,J)^SNiP). (2.12) 

Alternatively, one can consider also convergence of the quenched averages with respect the 
external noise 

E{fNW, J)) ^ fNiP), E{uNiP, J)) ^ «iv(/S), E{sn{^, J)) ^ sM. (2.13) 
For the sake of convenience, let us define the quenched average 

= N-^E{\ogZNiP,J)) = -pE{fNiP,J)). (2.14) 
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Clearly, a^iP) is convex and increasing in /3, for /3 > 0. Then we have the following 
elementary estimates (see for example [6]). 

Proposition 1. For the previously defined Q!jv(/S) the following bounds hold, uniformly 
in N, 

\og2<aNiP) <log2 + ^V4, for ^<;5 = 2Vlog2 = 1.665..., 

</3yiog2, for /3 > /3. (2.15) 

The proof (see for example [6]) is based on the following simple, but important, in- 
equality, relating quenched averages with annealed averages 

E{logZN{P,J)) < logE{ZN{f3,J)). (2.16) 

In the annealed average the J's participate to the thermodynamic equilibrium. A simple 
calculation shows 

<Tl...(TN 

= 2(]^) = (2exp^)^. (2.17) 

Therefore, for any ^ 

aiv(/3) <log2 + ^V4. (2.18) 
On the other hand, the lower bound in (2.15) is trivial, since convexity gives 

Z^(^,J)= exp(-^iy^((7,J))>2^exp(-^ ^ iy^((7, J)) = 2^. (2.19) 

(Ji---(Jn ai...aN 

As shown in [6], a simple argument, based on thermodynamic stability and the bound 
(2.18), gives the upper bound in (2.15) for P > In fact, introduce the function 

0(/5)=log2 + /5V4, (2.20) 

and define ^ such that, with (l)'{/3) = {dp(j)){^), 

m-Mm- (2.21) 

Therefore 

^ = 2Vfog2, 0(^) = 21og2, 0'(^) = v^. (2.22) 
Assume by absurd that there exists a point /3 > P where the upper bound is violated, i.e. 

aNiP) > ^Vfog2, 13 > p. (2.23) 
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Notice that by convexity of ckjv we would have at this point 



a'M > (^-^)-^(ajv(^) -«iv(^)) > {P-P)-\P^/^-2\og2) = (2.24) 

On the other hand, the average entropy is given by 

E{sn{I3, J)) = cxNiP) - l3a'^{/3), (2.25) 

as a consequence of (2.7,10,11,14). Therefore, by exploiting again convexity and the bounds 
(2.23) in P and (2.18) in P, we would get 

E{snW, J)) < c^nW) - PiP P)-^ (a7v(/3) - aNCP)) 
= -;5(/3 - B)~^aN{(i) + /3(/3 - ^)-^aN{^) 
< -P{/3 - ^)-l/3^/^^ + PiP - ^)-^2 log 2 = 0, (2.26) 

which is impossible. 

Of course, the bound given by Proposition 1 extends to the limit N ^ oo. 

Proposition 2. In the thermodynamic limit we have 

limsupajv(/5) <log2 + ^V4, for ^ < ^ = 20^ = 1.665 ... , 

JV->-oo 



<^Vlog2, for j3>/3. (2.27) 



As remarked in [6], these bounds show a phase transition for some /3c < /3 — 2^1og 2, in 
the sense that annealing of the noise J gives wrong results in the limit AT — >■ oo. On the 
other hand, next Theorem 3 shows that annealing is correct in the thermodynamic limit 
for any (3 < 1. Therefore, the transition must occur at some critical point (3c, such that 
1 < /5c < 2-\/log2. On physical grounds (see for example [2]), it is expected that /3c = 1- 
This is also confirmed by the buildig up of strong fiuctuations in the N ^ oo limit of 
global thermodynamic variables, as /5 t 1, shown in the beautyful and detailed analysis of 
Aizenman, Lebowitz and Ruelle [6], based on graph expansions. 
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3. THERMODYNAMIC LIMIT IN THE HIGH TEMPERATURE REGIME. 

Let us introduce the important order parameter M^{/3) (see for example [2], [6], [8]) 

mUp) = J_ . o<M|,(^)<i, (3.1) 

and write the average internal energy as 

-E{unW, J)) = c^M f (1 - MUP)). (3.2) 

This comes from a simple integration by parts on the external Gaussian noise, expressed 
in the form, reminiscent of the Wick theorem in quantum field theory, 

E{j,jF{J))=E{^^F{J)), (3.3) 

for any smooth function F of the noise. 

In fact, we can write, with obvious shorthand notations, 

-E{un{P,J)) = -N-^E{ujNiHN)) 

= {N^/N^)-^J2E{J,JUJN{a^aJ)) 

= pN-\N - l)-i J]£;(a;Ar(a,cT„cT,cT,)) 

= PN-\N - ^ (l - E{u;Ua,a,))). (3.4) 

We have exploited the general expression 

-^ujn{A) = {I3/^/N^)ujn{A, aiaj), (3.5) 

where un is the Boltzmann state ujn,i3,j, and oj{A,B) denotes truncation, i.e. oj{A,B) = 
u{AB) -u{A)u{B). 

Then we can state the following 

Theorem 3. With the definitions (2.14) and (3.1) we have 



lim aAr(/3) = log2 + /3V4, (3.6) 



lim Mf,{p) = 0, (3.7) 



N->-oo 



uniformly for < ^5 < 1. 



Here we give a streamlined proof, which exploits only elementary thermodynamic 
properties and positivity of fluctuations. We find convenient to split the proof in a long 
series of simple statements. It can be also understood as an elementary application of the 
cavity method. 

Let us define 

aN{f3,t) = N-^E{\og J2 cxp( Jija,a,)exp(-^^J,(7,)), (3.8) 

where we have introduced an additional parameter t,t>0, and an additional fresh noise Ji, 
i = 1, 2, . . . , A^, with the same properties of the noise Jij. Notice that aAr(/5, 0) = aN{l3), 
as defined in (2.14). Clearly, aN{P,t) is convex and increasing in P and t, separately. For 
the sake of simplicity, we call Ut the state associated to the new Boltzmannfaktor. When 
t = 0, ujt reduces to the even Boltzmann state lon. 
Let us introduce the order parameters 

M^(^,t) = iV-i^£;(a;2(ai)), 0<Ml{P,t) < 1, (3.9) 

i 

such that 

dpaM{P,t) = ^{l-M'Mt)), 

dtaM{f3,t) = ^{l-M%{l3,t)). (3.10) 

Notice that Mfj{/3,0) = M'^{/3) as in (3.1), while Mj^[/3) — 0, since ojn is even. In 
a sense, the parameter t gives a smooth interpolation between a system with N sites, at 
t = 0, and a system with N + 1 sites, at t = \/N/ {N — 1), with a small change in the 
temperature. The parameter t denotes the strength of the coupling between the original N 
spins and an additional (A^+l)th spin added, with Jj interpreted as Ji,7v+i, i = 1,2, . . . , N. 
A precise statement is given by 

Lemma 4. The following equalities hold 



MUPV{N-1)/N,P) = Mj,{/3^/{N-l)/N,p) = Ml,^M, 

{N+l)aN+i{P) = log2 + NaNiPViN-l)/N,p). (3.11) 

For the proof let us notice that we have M'^{p,t) = E(ojf{aia2)), Mj^{(3,t) = E[uj^{ai)), 
-^Ar+i(/5) = -E'(^Af+i(o"icTiv+i)) = E (u%_^i{aia2)) , due to the complete symmetry among 
the J's. In fact, for each of the alternatives in (...,...,...), we can write 

V (1,0-1, (71(72) exp(^= Y] Jij(7i(7j)exp(^= V J,ai) = 

ViV V-/V 
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(l,CriCr7V+l, (71(72) exp(^= JijUiUj), (3.12) 



2 ••• 2 ^ 

(Ji ...(JjVfTjv + i <Ti...crjv + i (*>i) 

where we have added the dummy summation variable ctn+i in the first equahty, and 
performed the change of variables C7j ctiCTn+i, i — 1,2, . . . , N , in the second equality. 

In the following theorem we exploit the mild dependence of aN{(3, t) on t, as a conse- 
quence of the 1/N factor in the t derivative in (3.10), and simple convexity estimates, in 
order to have bounds on M^{/3,t), uniform in t. 

Theorem 5. Assume /S > A/3 > 0, then we have, uniformly for < t < ^5, 

MKP + A« - M _ < M%(P, t) <MUP-m+f + (3.13) 

For the proof, let us notice that (3.10) gives 

< dta^iP, ^) ^ ^' < - aiv(/3) < (3.14) 

Therefore, by convexity we can write 

>aiv(/3)-aiv(/3-A/3)- 



2N 

>A/3a'^{P-A/3)-^. (3.15) 



^From this and (3.10) we have 



a'^{/3,t)>a'^{/3-A/3)- 



2NAP' 



^(1 - MUM) > ^^(1 - MUP - - (3.16) 



We conclude 



MUM < (1 - y-)MUP -Ap) + ^ + (3.17) 



This gives the upper bound in the theorem, and the lower bound is proven analogously. 
Through a straigthforward calculation and integration by parts, we can now prove 

2 

the following basic theorem, which gives the t derivative of the order parameter Mpj{f3,t) 
in terms of the order parameter M^{/3, t) and some fluctuations of variables of the theory 
(see also [8]). 
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Theorem 6. The following equality holds 

-E{{N-^Y.^l{a,)f)-AE{ut{A,A)), (3.18) 

i 

where A is the variable 

yl = Ar-i^a,a;i((7i). (3.19) 

i 

For the proof, let us evaluate 

at(iV-i J]£(a;2^a,))) = 2iV-i J] E(a;i(aOata;t(ai)) = 

i i 

= 2tN-^ ^ E(a;t((7i, aj) + ujtiai)ujtiai, aj, aj)), (3.20) 
where we have exploited 

1 d t 



UJticTi, (Tj) = —7=Ut{cri, (Tj, (Tj). (3.21) 



Now we consider the identities 

^^t{(^i, (^j, CTj) = -2a;t(cri, aj)ut{(Tj). (3.22) 
By collecting all terms, we can write (3.20) in the form 

dt^(N-'^^E{ul{ai))^ = J]£;(a;2^cTi(j,)a;2(cTi)a;2((j,) -M(cTi)a;t(cT,)a;^((7i,(7j)), 

(3.23) 

and the theorem is proven. 

Now we exploit the obvious positivity inequalities 

a;,(A,A)>0, i?((iV-i^a;2(a,))') >i?((iV-i^a;2(a,)))', (3.24) 

i i 
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and the inequality 



N 

= (1 - 1^)MUP, t) + j^<MU^,t) + ^, (3.25) 

together with the definitions (3.9), in order to derive 
Theorem 7. The following inequalities hold 

dt2Ml{^,t) + M%{p,t) < Mlip.t) + ^, for any t, 

< M%{I3, /3), uniformly for < t < /3, (3.26) 

where 

M''M~p) = max M|,(At) + ^, (3.27) 

M%{/3,t) < Mn{I3 J) tanh {t'^MN{/3j)), uniformly for < t < ^. (3.28) 

In fact, (3.26) follows from (3.18), (3.24), (3.25) and the definition (3.27), while (3.28) 
follows from a simple integration. 

Let us now put {N — 1)/N in place of j3 in (3.28), and then ^ — j3, and t — j3. We 

get 

MU/3V{N - 1)/N, p) = Mli^y/iN -1)/N, p) 

< MNiPViN - l)/N, 13) tanh {eMN{Py/{N - 1)/N, /3)). (3.29) 

Now we are ready for the proof of the following key theorem. 
Theorem 8. Define 

M^= max^M|,(/3), M|, < 1, (3.30) 



then for large N 



hence 

lim M|r = 0. (3.32) 

JV— >-oo 

For the proof, let us notice that the definitions (3.27) and (3.30), and the uniform 
upper bound (3.13) of Theorem 5 give for /3 < 1 



. , 1 .A^ViV , /5 , , 1 



<M^ + ^^=( + — ^) + — . (3.33) 
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It is convenient to put = P/y/N, so that 

MU^,/{N-1)/N,I3) < Ml, + -^^= + 1 (3.34) 
But in the same conditions we derive from the lower bound (3.13) 

MUPViN-l)/N,P) > M^(/3(V(iV-l)/iV + -^)) - ^i^. (3.35) 
Therefore, from (3.29), (3.34) and (3.35) we derive 

Mf,{-^{1 + VN^)) < -^^== + MVtanh(^2MV), 

Notice that (1 + ^/N — 1)/\/N > 1, therefore, if we take the maximum in (3.36), for 
< ^ < 1, we immediately have the following important bound 

M'l < ^j^=j + j^ + M'N tanh(^2^'Ar), (3.37) 

where we added 2/\/N — 1 + 1/A^ to both members for convenience. Since the function 
x{x — tanhx) is strictly increasing for a; > 0, and < M' n, we have also, 

Mjv(MAr - tanhMAr) < + ^, (3.38) 

which shows that Mjv as N ^ oo. 

In order to have an estimate on the rate of convergence, we exploit the following 
inequality, holding for < a; < 1, 

x{x - tanhx) > a;^(tanh 1)^/3, (3.39) 

and the theorem follows. 

Finally, we can prove Theorem 3. In fact, from M^{P) < Mjv, for < y5 < 1, we have 
the uniform convergence in (3.7). On the other hand, we also have 

«'iv(^) = f(l-M|,(^))>|(l-M|,), 

aNiP) > log2 + ^2(1 _ m2,)/4, (3.40) 

and the uniform convergence in (3.6) follows. 

In conclusion, we see that the validity of the annealed approximation holds uniformly 
in < /5 < 1, in the infinite volume limit, as a simple consequence of thermodynamic 
stability and positivity of the mean square fluctuations. 
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MARGINAL MARTINGALE AND CAVITY METHODS. 

Let a; be a generic even state on the Ising variables cxi, . . . , crjv, possibly depending on a 
stale noise Jij. Introduce the marginal martingale function iljN{u),t) defined by 

V'Tvl'^,^) = logo; (exp^= V^JjCTi), (4.1) 

V N 

where Ji are a fresh noise, as in (3.8). Since oj is even, we can substitute the cosh in place 
of the exp in (4.1). 

We have the bounds given by 
Theorem 9. For any t, the following bounds hold uniformly in uj 

J logcosh(tz) dfM{z) < ipN{i^,t) <N J logcosh{tz/VN) d^i{z), (4.2) 

where d^{z) = exp(— z^/2) dz/y/^m is the unit Gaussian distribution. The upper bound is 
realized if uo in (4.1) is taken as the symmetric product state where all configurations 
of the cr's have the same probability 2"-'^. This is the case if oj is the Boltzmann state ojn 
at ;S = 0. The lower bound is realized if oj in (4.1) is a state ID which gives equal weigths 
i to some fixed configurations ctl, . . . , cTjq and its inverted one — cti, . . . , — ctjv, and zero to 
all other configurations. 

Proof. First of all let us show that the bounds are in fact equalities for some states. 
For the upper bound let oj in (4.1) be the symmetric product state oj^^\ In this case we 
have 

loga;^°)(exp^ V JiCTi) = V log 0;^°) (exp ^J^cTi) = V log cosh(^ J^). (4.3) 
V-/V \J N V-/V 

By taking the average E and exploiting the symmetry in the J's, we have the upper bound 
in (4.2) as an equality. For the lower bound, let oj be the state oj which gives probability 
I to each of the two configurations CTj = 1, and = —1 for all z's, for example. In this 
case we have ^ ^ 

a;(exp —=Jiai) = cosh(^= Jj). (4.4) 

V A* V N 



By taking the E average we get the lower bound in (4.2) as an equality. 

Let us now prove the bounds for a generic state. Introduce the additional Ising 
variables ei, . . . , ejv, with symmetric distribution, and let E' denote the related average. 
Since the J's have symmetric distributions, by annealing in E', we have 

£;ioga;(exp ^= Jiai) = log a; (exp ^= Jtatei) 

< ElogujE'{exp ^= Jtaiti) 
VN Y 

= Vlogcosh(^J,), (4.5) 
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where we have freely exchanged E' and u). Therefore, we are reduced to (4.3), and the 
upper bound follows. Obviously, we could anneal completely the J variables and get the 
slightly weaker bound uniform in N 

£;ioga;(exp 5^ JiCTi) < (4.6) 

For the lower bound, it is enough to write 

£;ioga;(cosh -^^J^di) > a;(£; log cosh ^ J^di), (4.7) 



by quenching also the variables cr, in some sense. In the E average the variables a do not 
play any role, and the lower bound follows. 

Let us remark that the upper bound in (4.2) is realized for the state u^^^ of maximum 
entropy A?" log 2, while the lower bound corresponds to states U of minimum entropy log 2, 
compatible with eveness. It would be interesting to find the appropriate bounds if u> has 
some given fixed entropy. 

The interest in the introduction of the marginal martingale relies on the possibility to 
derive information on the thermodynamic variables from information on ipN{ui,t), for some 
properly chosen state. Let us remark that V'iv and related quantities play also a central 
role in the deep study of fluctuations made by Pastur and Scherbina [8,9]. We refer to 
[10] for an extensive treatment of the marginal martingale. Here we give only some typical 
results. 

Let us define the particular marginal martingale ipNiP), associated to the mean field 
spin glass model, by replacing the generic state u in (4.1) with the Boltzmann state a;^ 
at temperature /3^/ (N — l)/N, and putting t = 

^NiP) = i^N{uj'^^t) = N{aN+i{f3^l^^,/3) - aNW^I^^)). (4.8) 
^From (3.11) we have 



'N -1 

{N + l)aN+i{P) = log 2 + NaNiP^ — ^) + ^iv(/3). (4.9) 

We see that ipNiP) connects the free energy for a system with N + 1 particles with the 
free energy for a system with N particles, with a small change in the temperature. This 
is a typical feature of the cavity method (see for example [4]). 

The following theorem shows how information on the limiting behavior of iPn{/3) 
translates into information on the behavior of the thermodynamic variables. 

Theorem 10. Assume that a lower bound of the type 

Ml3)>^Nil3) (4.10) 
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holds uniformly ioi < (3 < P, N > K. Then we have 



liminfaAr(/5) >log2+ / ij{l3y/l^) dq, (4.11) 

for < /3 < p. An analogous statement holds for an upper bound. 
Assume that the following limit exists 

lim V'iv(/5)=V'(/5), (4.12) 

uniformly on a compact region < /5 < /5, with '0(/5) continuous in /3, as a consequence. 
Let us define 

«(^) = log2+ / dq 
Jo 

= log2 + r'/ i^{/3')dp'\ (4.13) 

^0 



so that the /3 derivative a'{/3) exists and the following holds 

a(/3)+/3a'(/3)/2 = log2 + V'(/3). (4.14) 
Then we have, ior < ^ < ^, 

lim aN{/3) = a{/3), lim a'j^{/3) = a'{P), lim ({N + l)ajv+i(/3) - NaN{f3)) = a{(3). 

N^oo N^oo AT— >-oo 

(4.15) 

Moreover, ipiP) ^tnd sixe increasing in (3, and 

W) = ^V2, (4.16) 

in the fluid region < /3 < 1. 

We give only a sketch of the proof, by referring to [10] for complete details. Let us 
iterate (4.9) N - K times 



{N + l)aN+i{/3) ^{N-K + l)log2 + KaK{/3\l ^^j^) 



N - 1 K 
+ ^n{(5) + ^N-i{(5\l + . . . + V'K(/3y -). (4.17) 



Therefore, by assuming the uniform bound (4.10), we have 



N 



K'=K 
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and derive (4.11) in the hmit N ^ oo, where the discrete variable K'/N becomes the 
continuous 1 — q. The other statements follow by standard arguments. 

A simple application can be obtained by combining the lower bound in (4.2) with 



This bound is a mild improvement over the analogous one found in [6] through a clever 
choice of trial states in the variational principle for the free energy. In fact, since log cosh a; 
is an even convex function of x, we have 



and the bound in [6] follows from (4.19). 

5. THE FUNCTIONAL ORDER PARAMETER. 

While the approximate solution of the model, given by Kirkpatrick and Sherrington in [1] , 
involves a numerical order parameter, function of the temperature, it was soon realized 
by Parisi that the complete solution [3,5] must be expressed through a functional order 
parameter x, which for a given temperature depends on an auxiliary variable q, both x and 
q taking values on the interval [0, 1]. This was suggested by numerical computations and 
a deep physical intuition about the rich structure of phases of the system, in the infinite 
volume limit, for /3 > 1. 

Here we prove that any marginal martingale can be expressed through a functional 
order parameter of Parisi type, also in the infinite volume limit N ^ oo. Moreover, we 
show that there is a very large set, in the convex space of functional order parameters, 
which gives equivalent results in the expression of the marginal martingale. This kind 
of gauge freedom in the choice of the functional order parameter raises the problem of 
picking the most convenient representation, according to some criterion. Parisi choice 
[2,3,5], fully discussed in [10] in the general frame introduced here, is very natural, because 
it gives simple physical interpretations for the parameters x and q, and very useful for the 
applications. But, in principle, other choices are also possible. 

The main result of this Section is given by the following 

Theorem 12. For a given even state w, on the Ising variables fxi, . . . ,(Tn, possibly de- 
pending on some stale noise, let us introduce, as in (4.1), the marginal martingale 





(4.19) 





(5.1) 



Then, for each oo and /3, there exists a functional order parameter 



X : 



[0,l]3q^x{q)e[0,l], 



(5.2) 
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such that 

W) = /(0,0), (5.3) 
where f{q,y), < q < 1, y & R, is the solution of the nonhnear antiparabohc equation 

(dJMy) + l{f"{q,y) + x{q)f'\q,y)) = 0, (5.4) 

with final condition 

f{l,y) = logcosh(/3|/). (5.5) 
In (5.4), /' = dyf and /" = ^yf- ^ possible expression of x is implicitely given by 

{l-x{q))E{u{f')) = N-'Yl^i^'(^^f'))^ (5.6) 

i 

where 

^ i ^ i 

and oj is the state 

u}{A) = uj{Aexpf)/uj{expf). (5.8) 

The proof is very simple, and can be understood as a consistent correction to successive 
annealing of the J variables [10]. It also shows how naturally the equation (5.4) arises from 
integration by parts. For a generic f{q,y) let us introduce the function defined by 

(j): [0,1] 9 g ^ (f){q) = E log w(exp /(g, V JiO-^)) 

VN . 

= Eloga;(exp/). (5.9) 

Notice the typical Brownian scaling ^/q in the y variable. 

Assume for / the boundary condition (5.5). Then we have 

0(l)=V'(/3), 0(O) = /(O,O). (5.10) 

Let us now calculate the derivative of (f){q). We have 

^ exp / = (dj) exp / + ^ ^ ^i^J' ^^P /• (^-H) 

Therefore 

Now we integrate by parts on the fresh noise Jj, by using 

^ exp / = .[la J' exp /, A/' = .flaj". (5.13) 
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By collecting all terms, we have 



m = maj) + 5/" + i/")) - ~ E^P'(<'*/'))- (5.14) 



I 



Therefore, if / is chosen so to satisfy (5.4) with x given by (5.6), we see immediately that 
(/){q) does not depend on q, and (5.3) follows from (5.10). Finally, let us notice that the 
inequalities 

0<uHaJ')<uif'^) (5.15) 

imply 

< x{q) < 1, (5.16) 

and the theorem is fully proven. 

We refer to [10] for a detailed discussion about the properties of the solution /(g, y; x) 
of the antiparabolic equation (5.4), with final condition (5.5), as a functional of a generic 
given X, as in (5.2). Here we only state the following 

Theorem 13. The function / is monotone in x, in the sense that x{q) < x(q), for all 
< Q < 1, implies f{q,y;x) < f{q,y;x), for any < q < I, y & R. Moreover / is 
continuous in the L^{dq) norm. In fact, for generic x, x, we have 

\f{Q,y;x) - f{q,y;x)\ < Y J \x{q') - x{q')\ dq' . (5.17) 

Let us now introduce the extremal order parameters xo{q) = and xi{q) = 1, such 
that for any x we have xo{q) < x < xi{q). It is simple to solve (5.4), (5.5) in these cases 
in the form 

MQ,y) = J logcosh{/3{y + zy/l^)) dfi{z), h{q,y) = logcosH/Sy) + /3\l - q)/2. 

(5.18) 

Therefore, we have for the / associated to a generic x the bounds 

fo{q,y)<fiq,y)<fi{q,y), (5.19) 
and in particular at the point q = 0, y = 0, 

/o(0,0) = J logcosh(/5;2) di^iz) < /(0,0) < {^^2) = /i(0,0). (5.20) 

We recognize that these bounds are the same as found in (4.2) and (4.6), written for t = /3. 
This remark leads to far reaching consequences. 

Theorem 14 Let ip{/3) be a marginal martingale for N particles, or a limit for A?" — >■ 00, 
therefore satisfying the bounds 



/ 



logcosh(^2) dii{z) < ipi^) < (5.21) 
18 



Let be a generic family of functional order parameters depending continuously in the 
norm on the variable e, < e < 1, with xq = 0, and Xi = 1, and nondecreasing in e. 
Then, there exists an e(/3) such that 

V^(/3) = /(0,0;x,(^)), (5.22) 

where / is defined by (5.4), (5.5) with x replaced by x^(^py 

The proof follows easily from the a priori bounds and continuity. An easy way to 
construct families of Xe is the following. Let x be some fixed functional order parameter. 
Introduce a^e, for < e < 1, defined by 

x,{q) = [x{q)-l + 2el (5.23) 

where [. . .] denotes truncation outside of the interval [0,1], i.e. [y] takes the value if 
y < 0, y if < y < 1, and 1 if y > 1. 
In conclusion, we have 

Theorem 15. If '(/'(/?) is a marginal martingale, then for each ^ there exists a nonempty 
hypersurface of functional order parameters such that 

^(/3) = /(0,0;a;), for any x G S^. (5.24) 
6. CONCLUSIONS AND OUTLOOK. 

We have seen that many important properties of the mean field spin glass model can be 
derived by elementary methods, based on the fluctuations with respect to the external 
noise. 

In particular, we have shown that the positivity of the mean square deviations alone 
forces the annealed approximation, for the free energy, to be correct in the thermodynamic 
limit, in the high temperature fluid region, up to the critical point. 

In general, the corrections to annealing can be easily expressed through a functional 
order parameter, of the type introduced by Parisi in his ingenious discoveries about the 
physical properties of the spin glass system. 

The methods exploited in this paper can be easily extended to other mean fleld disor- 
dered theories, with variables more general that the Ising variables, for example continuous 
ones. 

Finally, we refer to [10] for a more complete description of the martingale method and 
its applications. 
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